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Abstract 

We will examine the issue of diffeomorphism symmetry in simplicial models of (quantum) 
gravity, in particular for Regge calculus. We find that for a solution with curvature there 
$H ' do not exist exact gauge symmetries on the discrete level. Furthermore we derive a canon- 

ical formulation that exactly matches the dynamics and hence symmetries of the covariant 
picture. In this canonical formulation broken symmetries lead to the replacements of con- 
straints by so-called pseudo constraints. These considerations should be taken into account 
in attempts to connect spin foam models, based on the Regge action, with canonical loop 
quantum gravity, which aims at implementing proper constraints. 
\ We will argue that the long standing problem of finding a consistent constraint algebra 

for discretized gravity theories is equivalent to the problem of finding an action with exact 
diffeomorphism symmetries. Finally we will analyze different limits in which the pseudo 
constraints might turn into proper constraints. This could be helpful to infer alternative 
^\ [ discretization schemes in which the symmetries are not broken. 



1 Introduction 

In quantizing a given theory its symmetries play a crucial role. The question whether symmetries 
of the classical theory have also a representation in the quantum theory can have a drastic 
influence on the properties of the resulting quantum theory. 

For general relativity the symmetry in question is diffeomorphism invariance. As so far there 
is no satisfactory model of quantum gravity yet, also the fate of diffeomorphism invariance in a 
quantum theory of gravity is open. 

Nevertheless a successful implementation of diffeomorphism invariance into quantum gravity 
models could ensure the correct semi-classical limit and moreover help to resolve quantization 
ambiguities (see for instance pQ), that could otherwise render the models unpredictive. It is 
therefore important to discuss notions of diffeomorphism symmetries in the models at hand. 

A particular class of models, for instance Regge quantum calculus [21 [3], spin foam models [5J, 
(causal) dynamical triangulations [5], use discretizations of the underlying spacetime manifold 
as a regulator in order to define the specific model, i.e. a strategy how to perform the path 
integral. In particular in spin foam models the Regge action appears in a semi-classical limit 
[6]. We will therefore concentrate on the discussion of symmetries in Regge calculus. 

The main question regarding diffeomorphism invariance for such discretized models is whether 
a notion of exact diffeomorphism invariance can be found for the discrete model, or whether 
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exact diffeomorphism invariance can arise only in a continuum limit (alternatively in a sum over 
triangulations, that is all possible ways of discretizations). A notion of exact diffeomorphism 
invariance directly on the discrete level would simplify very much the process of defining the 
theory, for instance in choosing the path integral measure (which then has to be diffeomorphism 
invariant). Furthermore, as we will show in this paper, with such a notion it should be possible 
to find a canonical formulation for discrete gravity models, with a closed, i.e. consistent con- 
straint algebra. This is the main problem for canonical or Hamiltonian lattice gravity models 

mm- 

In this work we will show, that for models based on the Regge action (in 4d) , diffeomorphism 
symmetry is generically broken. This is contrary to expectations voiced in the literature [9]. 
Nevertheless there are many arguments that diffeomorphism invariance will be restored in the 
continuum limit pUl EH E2] , and also our results will support this view. 

These results do not exclude that discrete models with diffeomorphism symmetry can be con- 
structed. Indeed, the definition of symmetry that we apply, depends crucially on the dynamics 
of the model, as defined by the action, in this case the Regge action. Other actions might exist, 
that exhibit exact diffeomorphism symmetry. We will discuss such instances of different actions 
for the same system with exact and broken symmetries respectively for 3d Regge calculus with 
cosmological constant in section [6] and toy models in section [8l 

A long standing problem is the construction of a consistent discretized canonical model for 
gravity and a representation of diffeomorphism in such a model, see for instance |13|, I14j . In a 
canonical formalism gauge symmetries are reflected in constraints on the canonical data, that 
also serve as generators for these gauge transformations. As diffeomorphisms also include trans- 
formations of the time coordinate, general relativity is a so called totally constrained system. 
That is the Hamiltonian, the generator for the dynamics of the system, is a combination of 
constraints. Hence the constraint algebra is of central importance in order to have a consistent 
dynamics. 

Often canonical lattice models are defined by discretizing the constraints of the continuum 
theory. A typical problem in such discretized models is that the constraint algebra is not 
closed. This leads to severe problems for the quantization. According to the Dirac program 
constraints resulting from gauge symmetries have to be quantized and to be imposed onto 
the quantum states. That is however only possible if the constraint algebra is closed, or in 
other words anomaly free. Here the problem is already on the classical level, i.e. one has 
to face classical anomalies. Although different approaches exist to circumvent this problem 
[HU HU Q2] , it might be quite hard in these approaches to keep classical and quantum anomalies 
and ambiguities under control. This would be easier if we could construct discretized models 
with a closed algebra. Here we will show that this problem can be seen to be equivalent to 
finding an action, i.e. a covariant model with exact diffeomorphism invariance. If the action 
displays exact gauge invariance we will find constraints generating the gauge transformations in 
the canonical formalism. For actions with broken symmetries we will find a different picture: 
Instead of exact constraints - that is relations imposed by the dynamics that happen to involve 
data of only one time step, we will find pseudo constraints, that is dynamical relations which 
show a (weak) dependence also on the data on the next time step. This dependence can be 
interpreted as a dependence on the Lagrange multipliers lapse and shift. This allows to solve 
the pseudo constraints for lapse and shift. Hence gauge freedom is lost (as it is indeed broken on 
the covariant level) and there are no constraints left on the canonical data - a picture that was 
advertised in the consistent discretization program [18] . We will however argue that - as one 
expects gauge symmetries to be restored in the continuum limit - it might be more promising 
to keep the pseudo constraints. Indeed, data leading to solutions with a small discretization 
scale, are concentrated on a 'thickened constraint hyper surf ace'. In the quantum theory the 
connection between covariant models with exact gauge symmetry and canonical models with 
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constraints is, that the path integral acts as a projector onto the space of states satisfying the 
constraints [19]. For systems with slightly broken symmetries one would expect instead of an 
exact projector an approximate implementation of the constraints, similar to a delta function 
versus a Gaussian. 

The advantage of the technique used in this work is that the dynamics as defined by the 
covariant equations of motion and the canonical time evolution equations coincide and hence 
also display the same amount of gauge symmetry. These methods are in particular important 
for attempts to establish a closer connection between spin foam models and (canonical) loop 
quantum gravity. In the latter dynamics is based on exact constraints whereas in the former 
we expect diffeomorphism symmetry to be broken (as it is broken for the Regge action and the 
Regge action appears in a semi-classical analysis [6] of current spin foam models). With this 
different handling of symmetries by the two models the dynamics as defined by these models 
very likely also differs - at least on the discrete level. 

In the next section we will discuss the notion of gauge symmetries we are going to apply and 
how to test for the existence of these symmetries. Next we discuss an evolution scheme for 
Regge calculus, the so-called tent moves, which we are going to need for the construction of a 
solution with curvature as well as for developing a canonical formulation. We construct such 
solutions in section [5] and determine the eigenvalues of the Hessian for these solutions. In section 
E] we discuss a canonical formulation for 3d discretized gravity with cosmological constant using 
discrete actions with exact and with broken symmetries. Furthermore we analyze limits in which 
the broken symmetries might turn into exact ones. We will show in section [7]that if we have exact 
gauge symmetries in the action then we will find proper constraints on the discrete canonical 
data. Section [8] provides another class of simple examples of discretized theories with broken 
gauge symmetries. For these class of theories however one can always define a discretization 
with exact gauge symmetries. We will close with a discussion in section [H The appendices [A] 
and [B] contain a description of two physically different solutions with the same boundary data 
and formulas for geometrical quantities of simplices that we need in section [6] respectively. 

2 Definition of gauge symmetry 

In this section we will shortly describe the notion of diffeomorphism symmetry that we are going 
to apply. 

First of all we will consider continuously parametrized gauge symmetries. (As discussed 
in [20] there might also exist notions which are completely based on discrete transformations, 
such as a change of triangulation. We will comment shortly on a possible relation to contin- 
uum symmetries in the discussion section.) In the continuum such gauge symmetries lead to a 
continuous family of solutions to the equations of motions (with fixed boundary values for the 
variables), instead of just having one solution. We will apply the same definition to discrete 
models. That is we will speak of an (exact) gauge symmetry if the boundary value problem 
displays non-uniqueness of solutions, moreover this non-uniqueness should be parametrizable 
in a continuous way. 

As solutions to the equations of motions are extrema of the actions, this means that the 
action is constant in some directions exactly at these extrema. As already discussed in [20] it is 
important to check that these constant directions also persist at the extrema of the action, that 
is at solutions. (Away from solutions any direction perpendicular to the gradient is a constant 
direction.) This means that it is not sufficient to identify (possibly configuration dependent) 
transformations that leave the action invariant as these transformations might act trivially on 
solutions. In this case we will still have uniqueness of solutions. 

If there is a continuous parameter set of solutions then there exist directions at these solutions 
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in which the action is constant and also the first derivatives of the action are constant (and equal 
to zero). Hence in this case the Hessian, the matrix of second derivatives, of the action, will 
have null eigenvectors. This criterion is therefore a necessary condition for a gauge symmetry. 

Note that different solutions might have gauge orbits of different size. There might be 
theories with solutions with gauge symmetries and solutions without gauge symmetries. As we 
will see this is the case for 4d Regge calculus. There, flat solutions display gauge symmetries 
[2TJ [22l [23]. Vertices of the triangulation supporting this flat solution can be translated (in 
four directions) without changing the flatness of the solutions. In contrast, for solutions with 
curvature (and moreover without any flat vertices^) we do not find gauge symmetries. Hence 
we have a mixture of exact gauge symmetries (for flat vertices) and broken symmetries. 

We will show explicitly that the criterion of vanishing eigenvalues of the Hessian (evaluated 
at solutions) is violated for a 4d Regge solution with curvature. But this example will also 
show that gauge invariance is only slightly broken, hence we can speak of approximate gauge 
invariance. Namely we will see that some of the eigenvalues of the Hessian (per vertex) are very 
small as compared to the rest of the eigenvalues. Moreover in approaching the flat solution these 
eigenvalues go to zero (quadratically in the curvature). Hence analyzing the eigenvalues of the 
Hessian is a precise tool to discuss approximate gauge symmetries. This might be helpful in 
order to construct actions with an exact gauge invariance. 

In section [7] we will furthermore see that the criterion of null eigenvalues of the Hessian is 
related to the appearance (or non-appearance) of constraints in a canonical formulation of the 
theory. 

3 Tent moves 

Tent moves are a way of defining a discrete time evolution for Regge calculus, which has first 
been described in [24], and used for several works in Regge calculus [12] . It is a way to evolve a 
triangulated hypersurface locally, such that the triangulation (that is the adjacency relations) of 
the resulting new hypersurfaces does not change. These tent moves are a very convenient tool 
to define a canonical formalism [20] for Regge calculus. Implementing some ideas from discrete 
numerical integration [25] or 'consistent discretization' [18] the dynamics defined by the canonical 
and covariant formulations will exactly coincide and therefore also the gauge symmetries in these 
formulations. 

Consider a (d — 1) -dimensional triangulation S, which can be thought of as a triangulated 
Cauchy hypersurface. Pick a vertex v in the triangulation and define a new vertex v* lying in 
the 'future' of v, and connect both vertices with an edge. Denote all other vertices in £ that v is 
connected to by 1, . . . , n. Connect also v* to the 1, . . . , n by edges. Furthermore we will have a 
simplex v*ij(k) (with vertices v* , i, j in 3d and 4d and v* , i, j, k in 4d) in the evolved hypersurface 
£* for every simplex vij(k) in S. Hence the triangulations of the two Cauchy surfaces are the 
same. The evolution can be thought of as gluing a certain piece of d-dimensional triangulation 
onto the hypersurface. This d-dimensional triangulation consists of simplices vv*ij(k) for every 
simplex vij{k) in £ in addition to simplices in the boundary coinciding with either E, T,* or 
both. The edge connecting v and v* is called the "tent-pole" . 

Note that each tent move can be generated by a sequence of Pachner moves applied to the 
Cauchy surfaces. For d = 3, an n-valent tent move is the result of a 1 — 3-move, followed by 
(n — 2) 2 — 2 moves, and finally a 3 — 1 move. By applying tent moves to various vertices after 
another, one can build up a large d-dimensional triangulation. 

There are several advantages of this description: 

1 That is, all triangles adjacent to these vertices have vanishing deficit angles. 
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Figure 1: A 5-valent tent move at the vertex v. 



• The evolution is local, in the sense that only the triangulation containing the vertex v 
(called the "star of v n ) is evolved, while the rest of the triangulation remains untouched. 

• As a result, one can evolve a collection of vertices independently of each other if neither of 
the vertices can be connected to any other of the collection. This has in particular been 
implemented in numerical applications [12] . 

• The tent moves are particularly useful for an investigation of Regge calculus in a canonical 
language, since all Cauchy hypersurfaces that are produced in each step are isomorphic 
(as simplicial complexes). Therefore, in each step the number of canonical variables on 
the hypersurface remains unchanged. This is not true in an arbitrary triangulation, which 
makes the canonical analysis harder, and is the main reason why we consider tent moves 
in this work. Furthermore for the analysis of a tent move we need to consider only a small 
number of equations as opposed to a scheme in which all vertices are evolved at once. 

• The choice of which vertices of a hypersurface to evolve can be understood as a discrete 
choice of lapse (to be either vanishing or non-vanishing). Also, if v% and vi are vertices 
in the initial triangulation that are connected by an edge, the two tent moves applied to 
v\ and V2 do not commute. If one first "evolves" V\ and then V2, one obtains a different 
(ti-dimensional) triangulation than first evolving v% and then v\. That might serve as a 
starting point for a definition of a discrete notion for a hypersurface deformation algebra 

3.1 Evolution equations for tent moves 

The Regge action for a d-dimensional triangulation T with boundary dT and interior T° := 
T\dT is given [21 [26] by 

SRegge = - V h£h - J2v h ip h + \J2V* (3.1) 

heT° hedT aeT 

where the h are the d — 2 dimensional subsimplices (sometimes called "hinges" ) of T and a are 
the top-dimensional simplices of the triangulation. Vh and V a denote the volume of the hinge h 
and of the simplex a respectively. The deficit angles and exterior angles iph are given by 

e h = 2tt - °h for h ^ T ° ( 3 - 2 ) 
A = 7T - J2 e h f o r h G dT, 
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and in both cases the sum ranges over all d-dimensional simplices a which contain h, and 9^ is 
the interior dihedral angle in the simplex a between the two (d — l)-dimensional subsimplices 
that meet at h C a. 



Note that because of the boundary term in (|3.ip the action is additive if we glue two pieces 
of triangulations together. 

The equations of motion are obtained by varying the action (|3.ip with respect to the lengths 
l e of the edges in T°. The boundary term in (|3.ip ensures also that, for all edges e £ T° the 
Regge equations read 



EdVh . dV a 



hDe 







(3.3) 



The variation of the deficit angles appearing in the action vanishes because of the Schlafli identity, 
see appendix [Bj 

Consider a tent move for an n-valent vertex v in the boundary X of a d-dimensional trian- 
gulation T. After the tent move is performed we will have a new triangulation T* with (n + 1) 
new inner edges, namely the edges vi, i = l,...n adjacent in S to the vertex v and the tent 
pole vv* . We will denote by St, St* and S* the action (with boundary terms) of the original 
triangulation T, the new triangulation T* and the piece of triangulation added in the tent move, 
so that we have St* = St + S*. The equations for the new inner edges can then be written as 



dSr + dS 1 
dl v i dl v i 
dS* 








(3.4) 



With the definitions 



Pvi :z 

Pvv* *" 
Pv*i :z 



dSr 

dS T 

dl vv * 

dS, 

dl v *i 



the equations of motion (|3.4p are now given by 



P 



Pc 



Pi 



Pvi 
Pvv* 

Py*y** 

-'Pv 



dS* 
dl v i 
dS, 

_dS 1 







(3.5) 



(3.6) 



that is the momenta p~ and p + defined as derivatives of the actions associated to the two pieces 
of the triangulation T* have to coincide. Therefore we will omit the superindices +, — . With the 
second and third line in (|3,5p we use the action of the added piece as a generating function of 
first kind to define a canonical transformation from the canonical variables (l v {, lw* iPviiPw* 

) to 

a set of new canonical variables {lv*i^v*v** iPv*i^Pv*v**) (where we introduced a fiducial vertex 
v** which can be thought of as the vertex added in a second tent move). Here we see that 
the length of the tent pole has a special status as its conjugated momentum is constrained to 
vanish. Hence this variable is not fully dynamical. This corresponds to an analogous result in 
the continuum where in the canonical analysis the momenta conjugated to lapse and sift are 
constrained to vanish. 

Equations (13. 5|) are just a reformulation of the equations of motion in canonical language by 
using the action 5* as a generating function for a canonical transformation. This allows us to 
obtain a canonical formalism which reflects exactly the dynamics of the covariant formulation. 
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Other attempts to define a canonical framework for Regge calculus usually involve changing the 
dynamical set up [7]. Therefore gauge symmetries of the covariant formulation might not be 
reflected properly in the canonical formulation. 

The advantage in the formulation used here is that the dynamics defined in the canonical 
formalism is the same as the covariant dynamics. Gauge symmetries of certain or all solutions 
should therefore have repercussions for the canonical formulation. 

4 Remarks on discrete ambiguities of solutions to the Regge 
equations 

Given a boundary value problem the question arises whether the solutions are unique. Gauge 
symmetries in the form discussed here lead to a continuous family of non-unique solutions. In 
addition there might be discrete ambiguities. 

These also appear in Regge calculus, as for instance reported in [27J. In our investigations 
we noticed ambiguities already for the smallest boundary value problem, namely a triangulation 
with only one inner edge. Such ambiguities are common to discretizations of continuum theories. 
Typically there is only one solution that is useful in a continuum limit whereas the others can 
be seen as discretization artifacts. Nevertheless the question of discrete ambiguities should be 
explored in more detail, as these might influence the quantum theory. 

Some of the ambiguities arise because of the following: If we consider for instance a closed 4d 
ball with its 3d triangulated boundary, then pieces of this 3d triangulation might stick inwardly 
or outwardly. Similarly if we solve the tent pole equation (the second equation in (13. 4ft ) for 
the length of the tent pole then typically there is one solution which is forward pointing, i.e. 
a (bigger) tent is built on a smaller tent and another solution, that is backward pointing, i.e. 
resulting in the tips of two tents in opposite direction. Here we will always select the forward 
pointing direction for the canonical analysis. 

In the appendix Owe will discuss another kind of ambiguity (arising by setting the prefactor 
in the equations of motion (|3.3p to zero). There we construct a boundary value problem 
which allows for a flat and a curved solution (with the same boundary data). The curved 
solution has however very high curvature and one reason for its appearance seems to be the 
highly symmetric situation. It can therefore considered to be a discretization artifact. 

5 Construction of a Regge solution via tent moves 

Here we will describe shortly how to find numerically a small Regge solution with curvature 
using the tent moves. The evolution of four-valent vertices should lead to flat solutions^: indeed 
such solutions can be constructed by subdividing accordingly a flat 4-simplex. 

Therefore the simplest case to consider is the evolution of a five-valent vertex. To this end 
we have to define the three-dimensional triangulation around the vertex v, we want to evolve, 
more concretely the three-dimensional star of v . As v is five-valent we have five further vertices 
which we will denote by 1, . . . , 5. We will assume that we have six tetrahedra with vertices 

vl24, «134, v234, ul25, ,vl35, v235 . (5.1) 

Accordingly we will have nine triangles of the form vij with i,j = 1, ... 5 in this triangulation, 
five edges of the form vi and nine edges of the form ij (all possible ordered combinations of 
i, j £ {1, ... 5} with the exception 45). 

2 Apart from the discrete ambiguity described in the appendix. 
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Figure 2: The 3d boundary of the tent move triangulation. 



To simplify the calculations even further we will assume that all edges ij have the same 
length kj = s. As the Regge vacuum equations are invariant under a global rescaling we will set 
s = 1. We will also assume that l v { = a, i = 1, 2, 3 are equal to each other as well as l V 4 = = b. 
Hence we have to deal with two dynamical configuration variables a and b. By a n ,b n we will 
denote the values of these variables at the time step n. Together with the lengths of the tent 
poles t n these will be also the free variables in the action. (Varying the action with respect to all 
the inner edge length and then looking for solutions with l v ; L = a, l VK = b is equivalent to using 
this reduction in the action and varying with respect to a, b and the lengths of the tent poles. 
In this sense we will consider a symmetry reduced action.) 

The 4-simplices involved in the tent move are all of the same type v°v 1 ijn where u ,^ 1 
denote the two vertices of the tent pole (at time steps n = 0,n = 1 respectively), i,j take values 
in 1, 2, 3 and k in 4, 5. We will denote by 

9®, A® the dihedral angle and the area of the triangle v°ij, 

6®, A® the dihedral angle and the area of the triangle v°in, 

9f, A<1 the dihedral angle and the area of the triangle v°v 1 i, 

6%, A\ the dihedral angle and the area of the triangle v°v 1 k, 

Q\, A\ the dihedral angle and the area of the triangle v 1 ij, 

9\, A\ the dihedral angle and the area of the triangle v x ik respectively . 

The canonical equations of motion determining a 1 , 6 1 and t° given initial values a°,b° and 
p® and p® are given by 
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9A t ,„ ^ , , 9A t 



m0 (2vr - 46?) + 2^f (2. - 3^) (5.2) 

8A° f)A° f)A a 

Po = ^ (— ^ + (7r -2*g) + |£(2* -4%) 

P'o = 3M (7r _ 2^ + ^(2* -a*) . (5.3) 
The new momenta p\ and p\ are defined as 

A = -(^(x-!»8 + ^(»-M}) + ^(*r-«f)) 

A = -(#('- 2 "') + ^( 2 '- 39 ')) ' (M > 

To construct a solution with an inner vertex v l we have to consider at least two consecutive 
tent moves. We will proceed in the following way. First we will start with some initial data 
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(a , 6°, a 1 , b 1 ) and use equation (|5.2p to find the length of the tent pole t°. Given these five 



length we can determine the momenta (p\,pl) through equations (|5.4p . Now we have a full set 
of initial canonical data and can use all three equations (|5.2l5.3p to find t , a 2 , b 2 . We will end 
up with a Regge solution with inner vertex v . In the end we have to evaluate the matrix of 
second derivatives of the action with respect to (a , b 1 , t°, t 1 ) on this solution and determine its 
eigenvalues. 

However if one attempts to solve the equations (|5.2l5.3p for t 1 ,^,^ 2 numerically one will 
typically encounter the difficulty that standard numerical (iterative) procedures do not converge. 
This is due to the fact that there is at least an approximate gauge invariance (corresponding to 
choosing the lapse at the vertex). So if we reformulate this problem as finding the extrema of 
the action, we will have the difficulty that there exist a direction in which the action is almost 
constant and the extremum along this direction is difficult to locate. The corresponding small 
eigenvalue of the Hessian leads to convergence problems of the numerical procedures. 

This difficulty is usually circumvented by for instance (gauge) fixing the value for t 1 and 
ignoring one of the equations, for instance (|5.2p coming from the variation with respect to t 1 
|12j . One then hopes that this equation is satisfied up to a certain small error. This error is 
basically determined by the value of the smallest eigenvalue in the Hessian. 

To obtain a solution to a predetermined precision 5 we can guess a value for t , solve the 
second and third equations for a 2 and b 2 and then evaluate the right hand side of equation (|5.2p . 
If this error E(t\) is larger than 5 we have to start again with another value of t . This procedure 
can be systemized by starting with two values for t , say A, B and determining E(A), E(B). If 
the signs of E(A) and E{B) are different, than there is a zero in [A, B] (assuming that E is 
continuous on this interval). This zero can be determined iteratively to arbitrary precision. 

We want to generate a set of solutions deviating slightly from the flat solution. Hence we 
first construct a flat solution and then introduce a deviation from the flat boundary data. Flat 
solutions can be found by starting from a given data set (o°,6 ,^) and considering the deficit 
angles at the triangles v°v 1 i,i = 1,2,3 and v°v 1 k,k = 4,5. Setting these to zero gives two 
equations and we can find values a\(a°, b°, and frj-(a°, 6°, t®) satisfying these equations. Since 
we set certain edge lengths equal to each other we have in this case only one lapse degree of 
freedom, i.e. no shift degrees of freedom. This means that we find one flat solutions for any 
given value of t * (and a°,b°). 

Evaluating the Hesse matrix (second derivatives of the action with respect to 
(a 1 , b 1 , t , t 1 )) on flat solutions we will find one zero eigenvalue corresponding to the lapse degree 
of freedom. Hence for a solution with curvature we expect either one vanishing eigenvalue - 
in case that the symmetry of the flat solutions persists - or at least one very small but non- 
vanishing eigenvalue indicating that the symmetry is broken for curved solutions. 

To obtain curved solutions we start with boundary data (a°,b°,a^,bj + x) deviating from 
the flat solution. We solve numerically for t°, a 2 , b 2 and t 1 (with error terms of the order 10~ 13 ). 



-5.X10" 6 
-0.00001 

-0.000015 
-0.00002 

-0.000025 
-0.00003 



0.0002 0.0004 0.0006 0.0 



Figure 3: The lowest eigenvalue of the Hessian Figure 4: The deficit angle at the triangle v°v i 
as a function of the deviation parameter function of the deviation parameter x. 
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The plot in Figure is for initial data generated as described above with a = 6° = 1 and 
Oj = 1.16902, b l j = 1.13356 corresponding to the choice t* = 0.2. As can be seen from this plot, 
the lowest eigenvalue obtains non-zero values for solutions with curvature. Compared to the 
other eigenvalues the lowest eigenvalue is very small as can be seen for the set of data given in 
Table [TJ In Figure H] we plotted the deficit angle at one of the inner triangles as a function of 
the deviation parameter x showing that we are indeed considering solutions with non-vanishing 
curvature. Furthermore we see that the lowest eigenvalue seems to grow quadratically with the 
curvature. 



X 


eigenvalue 1 


eigenvalue 2 


eigenvalue 3 


eigenvalue 4 


0.0001 
0.001 
0.01 


1678.52 
1659.28 
1528.00 


570.127 
562.471 
509.495 


-394.796 
-383.002 
-306.545 


-2.91136xl0- y 
-2.87768 xlO" 5 
-2.57318xl0~ 3 



Table 1: Eigenvalues of the Hessian as a function of the deviation parameter x. The inner vertex 
has seven adjacent edges, but due to the symmetry of the problem, these are associated to only 
four variables. 

This shows that for these vertices with curvature we do not have symmetries which can 
be associated to translating the vertex. We constructed also other examples (see appendix El 
for one other solution) where again for solutions with curvature we do not find any vanishing 
eigenvalues of the Hessian. This does not exclude that some curved solutions exist, which have 
vanishing eigenvalues - however this seems to be a rather non-generic case. 

Apart from the result that the symmetries in Regge calculus are broken for vertices with 
curvature we want to point out that with the methods presented here we can easily check to 
which extent the symmetries are broken and how this scales with the curvature. The example 
considered here indicates that the symmetry breaking grows quadratically with the curvature. 
As we will see in section [6] this might allow for a canonical formalism in which we obtain proper 
constraints up to terms with a specific order in the curvature. 

6 Canonical analysis for 3d Regge calculus with cosmological 
constant 

In this section we will investigate the consequences of broken and exact gauge symmetries in 
the action for a canonical formalism. Here we will discuss 3d Regge calculus with and without 
a cosmological constant term as this example provides us with both cases, one in which the 
symmetries are exact and one in which they are broken. Moreover it allows us to discuss limits 
in which the broken symmetries may become exact. The analysis of 4d Regge calculus (in which 
the symmetries are generically broken) will be postponed to a seperate paper [28J. 

The solutions of 3d dimensional general relativity without a cosmological constant are locally 
flat. Correspondingly the 3d Regge equations just require that the deficit angles, which are 
attached to the edges, should be vanishing. Hence every triangulation of a locally flat space is a 
solution and moreover there is a three-parameter gauge freedom attached to every inner vertex 
of the triangulation. This gauge freedom corresponds to translating an inner vertex such that 
the triangulation remains flat. 

This gauge freedom vanishes if we consider a non-vanishing cosmological constant and use 
the standard Regge action with an added volume term 

eeT° eedT t£T 
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where l e , e e are the length and the deficit angles respectively associated to the edge e, tp e is the 
exterior angle at a boundary edge e G dT and V^- is the volume of a tetrahedron r [291 120j . 

Indeed we can consider also in this case small 3d triangulations T with an inner vertex 
and evaluate the Hessian on a solution. In this case we used two consecutive tent moves at a 
three-valent vertex to construct such a small 3d triangulation. The results show small but non- 
vanishing eigenvalues corresponding to the approximate gauge symmetries, see Figures [5] and [6) 
We considered homogeneous solutions for which the three length variables at the evolved vertex 
coincide. The initial values are 1° = 1, 1 = 1.8 for the length variables at time step and 1, and 
s = 1 for the non-dynamical edges not adjacent to either v , v , v 2 . Nevertheless we can consider 
the Hessian with derivatives for all five variables (the three length variables and the two tent pole 
variables) and evaluate this Hessian on the homogeneous solution. Because of the homogeneous 
configuration there are two small eigenvalues coinciding which correspond to translations of 
the vertex in space- like directions (associated to spatial diffeomorphism constraints), that is 
variations for which the tent pole variables are constant. Then there is another small eigenvalue 
corresponding to a translation in tent pole direction (associated to a Hamiltonian constraint). 
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Figure 5: The eigenvalue corresponding to 
translations in tent-pole direction function 
of A = 0.01 x x-axis label 



Figure 6: The eigenvalues corresponding to 
translations in space-like direction function 
of A = 0.01 x x-axis label 



What makes the 3d case with cosmological constant interesting for us is that there is actually 
an action, defining a discretized dynamics, in which these gauge symmetries are exact [30) . 
Starting with the standard 3d Regge action we can therefore test different methods to construct 
an action with exact gauge symmetries and see if we obtain the same result. 

In this section we will start with the 3d standard Regge action and perform a canonical 
analysis using tent moves. Since the gauge symmetries are only approximate we will find that 
we do not obtain exact constraints. Rather we have to deal with pseudo constraints which 
depend on lapse and shift. 

One could try to obtain exact constraints by considering a limit in which the length of the 
tent pole (and therefore lapse and shift) goes to zero. Indeed such a strategy works for discretized 
reparametrization invariant systems as is explained in section [8) In this case however one cannot 
obtain useful results. The reason is that gauge invariance is not being restored in this limit. 
Note that here we keep the spatial discretization scale fixed and send the discretization scale in 
time direction to zero. That results in almost degenerate simplices. As is also argued in [llj one 
cannot hope for a restoration of diffeomorphism symmetry for such degenerate simplices. 

On the other hand one can consider the limit in which both the discretization scale in spatial 
and time directions are small. By rescaling the edge length and the cosmological constant 
appropriately this corresponds to small cosmological constant A. We will see that to first order 
in A the constraints do not depend on lapse and shift. Moreover these constraints correspond 
to the first order constraints obtained from the alternative action with exact gauge symmetries. 
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6.1 Canonical analysis of a tent move 



We will consider the tent move evolution of a three-valent vertex v. We will denote the length 
of the edges between v n , where the superindex n denotes the time step, and the adjacent 
vertices 1,2,3 by If, i = 1,2,3. The length of the tent pole between vertices v n and v n+1 is 
t n . Furthermore the (non-dynamical) lengths of the edges ij will be denoted by su. With 
S( n ,n+i) we denote the action with boundary terms for the piece of triangulation that is added 
by performing a tent move. In this case the added piece are three tetrahedra joined at the tent 
pole. 




Figure 7: A tent move for a three-valent vertex. 
According to (|3.5l3.6p the canonical time evolution equations are given by 

„ _ dS {n,n+l) = ,„+ _ dV(n,n+l) 
Pi Ql n Vi A dl n 

n _ 5S^+i) _ 9V(n,n+l) 
Pt ~ Q t n ~H Q t n 

n+ l ^(n,n+l) _ („+!)- dV {n , n+l) 

Pi m n+l % + A dl n + l 

Pt +1 = 9 -%£^ = ■ (6.2) 

Here tp™ + is the exterior angle at the edge v n i of the triangulation between time steps n and 
(n + 1) whereas we denoted with ^ n+l ^~ the exterior angle at the edge v n+1 i of the same piece 
of triangulation. is the deficit angle at the tent pole v n v n+1 and Vi n , n +\) the volume of the 
triangulation between the time steps n and (n + 1)- 

Because of the last equation in (|6.2|) we have that pf = for all time steps n. We want to 
solve the second equation p™ = for t n = t n (l™, l™ + ) and use this solution in the first equation. 
For A this equation however involves sums of trigonometric functions and cannot be solved 
explicitly. Nevertheless it is possible to make an ansatz t n = °t n + A l t n + . . . and solve this 
equation order by order in A. 

The zeroth order gives the dynamics for vanishing cosmological constant, that is solutions 
are locally flat and all deficit angles have to vanish. The solution °t n to the second equation 
can be easily constructed by geometrical considerations, it is the length of the edge between the 
two tips of a double pyramid with triangular base. (In general there are two solutions, one in 
which the tips of the pyramid are in opposite direction and one in which the tips are in the same 
direction. We will use the latter case.) Using this solution in the first equation in (|6.2I) (with 
A = 0) one will find that the dependency on the length drops out and that the momenta 
are constrained to be 

p? = -ir + 9 vni (l?, Sij ) (6.3) 
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where with 9 v ni(lf, sy) we denote the dihedral angle at the edge v n i of a tetrahedron with edge 
lengths If and Sij and with vertices t> n ,l,2,3.. Indeed since the dynamics of the theory just 
results in flat space the intrinsic and extrinsic geometry of the hypersurface resulting from the 
tent move evolution will be just the same as that of the boundary of a flat tetrahedron embedded 
in flat 3d space. The constraints f|6.3|) just express this geometrical fact. 

If we instead use the solution °t n = °t n (lf, lf + ) in the third equation of (|6.2[) . which defines 
the momenta pf +1 at the time step (n + 1) we will find that the dependence on the length 
variables If drops out and that the constraints are preserved by time evolution^] 

p -+l = _ 7r + ^ n+1 .(in+l, Sii ) . (6.4) 

To summarize, for A = 0, the canonical evolution equations (|6.2|) do not determine uniquely 
the evolution for the lengths lf +l and the momenta pf +1 , they rather result in three constraints 
determining the momenta pf and pf +l as function of the lengths If and lf +1 respectively. 
Therefore given three length If the momenta pf are determined by the constraints, whereas the 
three length lf +l at the next time step can be chosen freely (respecting generalized triangle 
inequalities). This corresponds to a free choice of lapse and shift at the vertex v n . The set of 
six length variables If, lf +1 determines the length of the tent pole t n . 

Moreover one can check explicitly that the constraints f)6.3[) are first class, even Abelian (see 
appendix [B]) . See also [31] , where it is shown that taking into account constraints based at 
different vertices we obtain a first class algebra. The constraints are the infinitesimal generators 
for time evolution - the free choice of the three Lagrange multipliers corresponds to the gauge 
choice of the three edge lengths in the equations (|6.2p . The choice of only one of the three 
Lagrange multipliers Nj = e for the constraints C, = pf + it — v n i(lf, s«) to be non-vanishing 
corresponds to the choice Zjwj = 1%^, = lj + £ for the lengths at the next (infinitesimal) 
time step. 

If we switch to a non-vanishing cosmological constant the solutions of the Regge equations 
do not display the same gauge freedom. This has consequences for the canonical analysis. 

The difference to the flat case is that if we solve (numerically) the second equation of (|6.2|) for 
the length of the tent pole and use this solution t n (lf, lf +l ) in the first equation the dependence 
of the momenta pf on the lengths lf +l will not drop out. Hence we can solve these equations 
for the lengths lf +l at the next time step as a function of the lengths If and momenta pf at 
time step n. Using the last equation in (|6.2[) we can then determine also the momenta pf +1 at 
the next time step (n + 1). That is given some initial data l®,pf the evolution is unique for all 
other time steps (ignoring discrete non-uniqueness and assuming the existence of a solution). 

Note however that the dependence of the momenta pfil™,^^ 1 ) on Z" +1 is very weak (for 
small A), see Figure El Indeed as we will see later, if we expand pf(Vj, in A only the second 
and higher order terms will depend on ■ 

Figure [8] shows the functions pf(lj',l 1 j +l ) evaluated on a homogeneous configuration where 
I'j = 1° and l r - +1 = I 1 as well as Sy = 1. We plotted the momenta as functions of 1°, but for 
different values of I 1 ranging from 1° + 0.01 to 1° + 1.2. This is a wide range considering that Iq 
varies between 0.8 and 2. Moreover we used A = 0.6, 1.2 and A = 2.4. The four plots at the top 
in Figure [8] are for the highest value of the cosmological constant A and show quite a variation 
with I . For A = 1.2 the variation is much smaller (the family of plots below the four top plots). 
The plots for I 1 = I 2 + 0.01 and I 1 = 1° + 0.4 are almost above each other. Choosing canonical 
data outside of the region traced out by the plots would lead to an evolution with large edge 

3 This can be seen from the fact that S(„ in+ i) is symmetric under exchange of and The different signs 

in the definition of momenta in (|6.2p are absorbed by the property of the momenta to change sign for evolution 
in 'forward' or 'backward' direction. See also the discussion for the momenta to first order in A. 
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Figure 8: The momentum p° as a function of the length 1° for A = 0.6,1.2,2.4 and I 1 — 1° + 
0.01,0.4,0.8,1.2. 

length or would not allow for any (real) solutions. Finally for A = 0.6 all the four plots are 
almost above each other. 

That is although we can apparently freely choose initial data l®-,p® a corresponding solution 
might not exist or lead to very large edge lengths. If one wants for instance the lengths l\ at the 
following time step to be in rather small intervals around (such that the length t° is small), 
the momenta will be also quite restricted to a small interval. For the flat case we argued that 
the gauge choice of the corresponds to a choice of the Lagrange multipliers lapse and shift. 
If we consider a non-vanishing cosmological constant we can interpret the functions 



as lapse and shift dependent constraints. For lapse and shift in small intervals we will obtain 
instead of one constraint hypersurface of infintesimal width a family of hypersurfaces labeled by 
lapse and shift or in other words a thickened constraint hypersurface of finite width. In the limit 
of vanishing cosmological constant (but keeping the lapse and shift intervals fixed) the width of 
this thickened constraint hypersurface converges to zero. 

Next we will discuss limits in which one may regain proper constraints. To obtain such 
constraints for finite A one might try to get rid of the l" +1 -dependence of the pseudo constraint 
Cf in (|6.5p by considering the limit Z™ +1 — > If (in case it is well defined). It corresponds to a limit 
in which the discretization scale in time direction goes to zero (whereas the discretization scale in 
spatial direction is fixed). Such a limit does indeed work for discretizations of reparametrization 
invariant systems and moreover is suggested by the plots in Figure El 

For the system considered here these considerations do not lead to useful results however. 
A reason for this can be found by considering the behavior of the eigenvalues of the Hessian of 
the action. That is we consider two consecutive time steps and consider the Hessian associated 
to the inner vertex. In Figures I9|10l we plot the eigenvalues of this Hessian, that correspond to 
translation of the inner vertex in time like and space-like directions respectively as a function 
of the length of the second tent pole. (The cosmological constant is A = 0.5, and we again 
consider a homogeneous solution with I® = 1, = 1.) The eigenvalue for the time- like direction 
- corresponding to the Hamiltonian constraint - does indeed go to zero with the length of the 
tent pole. The eigenvalues corresponding to the spatial diffeomorphism constraints however start 
rather to grow for very small lengths of the tent pole. Hence in this kind of limit only one of 
the three gauge symmetries seems to get restored □ This is in agreement with other arguments 
|llj claiming that for a restoration of diffeomorphism invariance for Regge calculus one needs to 

4 Note that if we consider a "symmetry reduction" of the system, in which all the Sij are the same as well as the 




= v n i-v\ 



(6.5) 
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consider triangulations in which the simplices are "fat" enough, i.e. are not degenerated. But 
in the limit we are considering here the tetrahedra become infinitely thin in time- like directions. 



0.006 - 
0.005 - 
0.004 - 
0.003 - 
0.002 - 
0.001 - 



Figure 9: The eigenvalue corresponding to Figure 10: The eigenvalues corresponding to 
translations in time-like direction as function of translations in space-like direction as a function 
the tent pole variable. of the tent pole variable. 

For the dynamics defined by the standard 3d Regge action with a non-vanishing cosmological 
constant we do not obtain constraints by trying to take the limit of a vanishing discretization 
scale in time like direction while the spatial discretization scale is kept fixed. (For actions with 
an exact gauge invariance such a limit can be trivially obtained as there arbitrary values for the 
lapse functions are allowed and hence also infinitesimal values.) We can also consider the limit 
in which both the time-like and space-like discretization scale is very small. By rescaling the 
length variables it is easy to see that this corresponds to taking the cosmological constant A to 
be small. Hence to consider this limit we can perform an expansion of the constraints in A. This 
expansion also allows us to find the constraints to a certain order explicitly. 

The zeroth order of the constraints is given by the flat dynamics (|6.3p . To find the first order 
we will make the ansatz t n = °t n + A x t n in the second equation of (|6.2h . expand everything to 
first order and solve for t n . We use this result 

4"=©"^ (6.6) 
\dt n J dt n \t"=H" 

in the first equation of (|6.2|) and expand again to first order in A: 

Pi __ 7r + < WZ ., Sv) + A 1^1 I +0( A ). (6.7) 

Although the explicit expression (|6.7p as function of the length variables looks quite lengthy (see 
appendix [B] detailing the derivatives of simplex volumes and dihedral angles) the dependence of 
the first order term on the length variables Z™ +1 at the upper time step drops out. This can be 

seen by rewriting the derivative —q^t- in (|6.7p as 

v"v n+1 ij p,nv n v n+1 ik 



v n v n+l ik 



V' L V 



del 



(6.8) 



I" we reduce the (potential) gauge freedom to translations in time direction. This is the situation used for the plot 
in Figure [S] In this case one can indeed find the time continuous limit and hence a time reparametrization invariant 
system, as the gauge freedom for these translations is restored for infinitesimal time steps. However beginning with 
the second order in A the resulting Hamiltonian does not coincide with the corresponding Hamiltonian constraint 
obtained from a symmetry reduction of the alternative dynamic with exact gauge invariance, see below. 
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where we used that = for the dihedral angles 01, 8^, at the edges e, e' respectively of a 
tetrahedron r (see appendix [B|) . In (|6.8p the tetrahedron r is the one with vertices v n v n+1 ij or 
v n v n+1 ik and e the edge between vertices v n i and e' between v n v n+1 . 

Now if we denote by °t the function of the variables 1™ , Z™ +1 that results by solving the 
equation e™(t n , Z™, Z™ +1 ) = for t n , we find by taking the total derivative of this equation with 
respect to If that 



(6.9) 



d°t _ def (defy 1 
dlf = ~dlf \d^J |tn=0t , 

Using (|6.8p and (|6.8p we have for the first order term in ([67 

/ (defy 1 gV^n+1) aV(n,n+l) \ _ / 

[dtn {dt"J at- " azf ; |tn=otn V ^ 9z™ az? , |tn=0tn 



dv (n , n+1) (°t(i?,q +1 ),q,i 



dif 

(6.10) 

The volume V (n)n+1) (9t(Z?, Z™ +1 ), Z",Z" +1 ) is the one of a flat "double pyramid", depicted 
in Figure [7J For the kind of orientation of the two tetrahedra as shown in this figure this 
volume is just the difference of the volumes of the two tetrahedra with length Z™, and Z™ +1 , Sjj 
respectively. Therefore the final expression for the momenta to first order in A is then 

dV(l n s- ) 

pf = -7r + 6 v n t (lj, Hi) + A y p %3) + 0(A 2 ) (6.11) 

where V(lj",Sij) is the volume of a tetrahedron r with edge lengths (Z",Sy). Note that the 
constraints truncated to first order are also first class and even Abelianjfl This can be easily 
seen by realizing that S T = Yleer^i^ — #e) ~~ AI/ T 15 a generating function for the first order 
momenta, that is pf = —dS T /dlf + 0(A 2 ). 

Furthermore the constraints are preserved under time evolution. Note that the action 
S( n ,n+l) is symmetric under the exchange of the variables If and Z™ +1 . For the first order 
of the momenta pf +1 as defined in the third equation of (|6.2p we therefore obtain 



1 n+l _ *V{n,n+l)(°t(lj,lj + „ , 



i 



that is again the derivative of the volume V( n n+1 ) = V(l™ + ,Sy) — V{lj,Sij) of a flat double 
pyramid (but this time with a plus sign). We have however to take the derivative with respect 
to the length lf +1 , that affects the larger tetrahedron of this double pyramid. In the end we 
obtain the same sign as in (|6.1ip 



pf^ = - n + 9 vn+li{ lf+\ SlJ ) + X --J— r ^ + 0(X') . (6.13) 



dV(q +1 



The constraints (|6.12p coincide to first order with a first order expansion in A of the con- 
straints describing an exact discretization of 3d gravity with cosmological constant [30]. The 



5 Note that we did not check the commutator between constraints based at neighbouring vertices. If the 'Cauchy 
surface' we are considering is the surface of a tetrahedron, we would get the same first order constraints for the 
four three-valent vertices and these would still be Abelian. 
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continuum solutions of this theory are spaces with homogeneous and constant curvature deter- 
mined by the cosmological constant. Accordingly instead of a flat tetrahedron embedded in flat 
space, we consider a tetrahedron with homogeneously curved geometry embedded in a space with 
the same kind of geometry. In analogy with the constraints (|6.4p for the flat geometry, which 
fix the momenta to agree with the dihedral angles of a flat tetrahedron, the constraints in this 
case fix the momenta to agree with the dihedral angles of a homogeneously curved tetrahedron 

OTad p? = -7r + A «U(*?,%) (6.14) 

(see appendix[B]for expressions giving the dihedral angles x 9 e for a tetrahedron in homogeneously 
curved space as a function of the edge lengths) . An expansion of these constraints to first order 
in the curvature or cosmological constant gives (|6.1ip . To see this one needs the identity 

dA " dl e ' ( j 

where X V is the volume of a tetrahedron in homogeneously curved space, that will be derived 
in appendix iBl 

Starting with second order in A the (pseudo) constraints (|6.5|) do depend on the length 
variables Z™ +1 at the upper time step. (This can be more easily checked by considering the 
"symmetry reduced" theory, where all Sij = s and all = l n , see also [20 1.) 

To summarize, we have seen that the canonical equations of motions for the tent moves 
reflects the gauge symmetries of the covariant theory. If the symmetries are exact, we will 
encounter proper first class constraints. (See also next section for a general derivation.) These 
constraints generate translations of the evolved vertex and hence mirror the non-uniqueness 
of the covariant solutions. If the gauge symmetries are broken, we do not obtain constraints 
in the usual sense, as the expressions depend on the lenght variables on the next time step, 
or equivalently on lapse and shift. This again mirrors that the covariant solutions are unique 
(ignoring discrete cases of non-uniqueness) and hence lapse and shift are fixed by the pseudo 
constraints. However to obtain reasonable (small) values for lapse and shift, the initial data have 
to be chosen carefully - effectively from a "thickened" constraint hypersurface of finite width. 

Trying to construct a time continuum limit in order to obtain proper constraints however 
fails in the case of 3d Regge calculus with cosmological constant. The reason is that not all 
gauge symmetries are restored in this limit, in which the time discretization scale goes to zero 
whereas the spatial discretization scale is fixed. 

If we take all edge length to be small, or equivalently consider a small cosmological constant, 
we can perform a perturbation in the cosmological constant A. The dynamics truncated to 
first order has gauge symmetries and we can obtain constraints. For 3d Regge calculus with 
cosmological constant exact constraints (reflecting exactly the continuum dynamics) exist [30] 
and are a higher order continuation of the first order constraints derived here. 

We will show in [28J that a similar expansion in curvature is also possible for 4d Regge 
calculus. In this case we do not know the exact (discretized) constraints but if we are able to 
construct the first order constraints the question arises if higher order terms can be derived and 
whether these are determined uniquely. These questions will be subjects for further research. 

7 Relation between symmetries of the action and constraints in 
the canonical framework 

Here we will discuss the relation between gauge symmetries of the action and constraints in 
more detail. 
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We argued that gauge symmetries should lead to a non-trivial action of the gauge group on 
solutions. As we are considering continuous groups this should result in non-unique solutions 
parametrized by the gauge group parameters. Since solutions are equivalent with extrema of the 
action, rather than having one isolated extremum there is a submanifold of extrema on which 
the action is constant. As the first derivatives of the action vanishes by definition we obtain as a 
necessary condition for continuous gauge symmetries that the Hessian of the action should have 
null directions. 

We will consider a triangulation obtained from two consecutive tent moves. Hence we will 
have edges with length Z™ -1 and in the 'lower' and 'upper' boundary of the triangulation 
respectively. Here e is an index for the edges from the evolved vertex v to the adjacent vertices. 
At the inner vertex v n will hinge edges with the length Z™ , i n_1 and t n . Denoting with S = 
S( n -i,n) + ^(n.n+i) the action with boundary terms for this triangulation the requirement for a 
null vector Y a , a = e n , t n ~ l ,t n gives the equations 

Eye ® 2 S , yt 71 - 1 9 2 S t n d 2 S 



dl%d1% dt n ~ x dl d dt n dl e ' 

\ ye 1 yt 



d\V:dt n - 1 dt n ~ l dt n ~ l 

e e 

YY e 92 S +Y tn 92 S = . (7.1) 
^ dl n e dt n dt n dt n y ! 

The momenta conjugated to the tent pole variables have to vanish and the length of the tent 
poles do not appear as boundary data. Hence these variables are not fully dynamical. Therefore 
we will first integrate out these variables. 

We solve for the lengths of the tent poles i™ -1 and t n the equations of motions 



dt n - 1 

^>(n,n+l) 





. (7.2) 



dt n 

Using the solutions r n_1 (Z^" 1 , l n ) and T n (Z",Z™ +1 ) in (j?T2"j) and differentiating these identities 
with respect to Z™ we obtain 

d 2 S {n _^ n) d 2 S (n _ hn) dT™- 1 
dt n ~ l dl™ dt^dt™- 1 <9Z™ 
d 2 S {n , n+l) d 2 S {n , n+l) dT n 
dt n dl™ dt n dt n dl™ ' 1 ' ' 

The resulting effective action is 

S{l n e -\r e J n e +1 ) :=S{l n e -\l n el l n e + \T n -\r e -\l n e ),T n {l n e ,l n e +1 )) • (7.4) 
Using the equations (17. 3p the Hessian H of the modified action S can be written as 
d 2 S d 2 S | d 2 S dT n - x d 2 S dT n 



dl?dl" dttdl n . dttdt n - x dl n . dVidt n dl n . 



Finally with the help of the three equations in (17. ip and again f)T.3f) it follows that null vectors 
for the Hessian H of S define also null vectors for the Hessian H of S: 

s^Hr ■ < 7 - 6 > 
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Because of the second and third equations in (|7.ip we can express the components Y tn 1 and Y*™ 
as a combination of the components Y e . (Here we assume that the second partial derivatives of 
the action with respect to i™" 1 and t n do not vanish. This is generically the case for the Regge 
action.) For this reason a set of linearly independent null vectors for the Hessian H will define 
a set of independent null vectors for the Hessian H of the same size. 
We will use the action S to find stationary points 

mr\ie +1 ) (7.7) 

as functions of the boundary data ~ , l™ + . Again we use these solutions in the equations of 
motion 

^L(%-\L n e (ir l ,l n e +l )l n+1 ) = (7.8) 
and differentiate these identities with respect to C, : 

^ d*S dl% Q 2 %, w+ 1) _ , . 

y di™di n e , di n + l di^di^ 1 ' { ' } 

That is we have 

H-L = -K (7.10) 



where L ee , := dL^/dl^ 1 and 



di?dr, +1 v ' 



Therefore null eigenvectors for the Hessian H are left null eigenvectors for the matrix K. 

Note that these considerations show that K has at least as many null vectors as the Hessian 
H . It is not excluded that K has even more null vectors, however that was not the case in the 
examples we investigated. 

The matrix K appears as the matrix of derivatives of the momenta 

Vl := -^§p (7.12) 

with respect to the length variables ■ If there exist null vectors we cannot expect to be 
able to solve the equations (|7.12p for the lengths +1 at the next time step as a function of the 
momenta and lengths at time n as the conditions of the implicit function theorem are violated. 

If we assume that the rank of the matrix K is constant (and not maximal) in a neighbourhood 
of the initial data set under consideration then we can conclude that the Legendre 

transform from the direct product of configuration spaces Q n x Q n+l to the cotangent space 
T*(Q n ) defined by 

(M +1 )~(M) (7-13) 

is neither injective nor surjective. 

The image of this map into the phase space is the (primary) constraint hypersurface. The 
number of (irreducible) constraints describing this surface is larger or equal to the number of 
linearly independent null vectors for the Hessian of the action in (|7.ip . The constraints are 
relations between the momenta and length variables at the same time step n and will follow 
from the equations (I7.12|) . 

This shows that if the action has (local) gauge symmetries, we will obtain proper constraints 
on the canonical data. (For general actions there might also appear constraints without having 
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gauge symmetries.) Also we cannot solve the relation (|7.12p uniquely for the length as 
a function of the data l^,Pe a t time step n. Hence there will appear arbitrary parameters 
(Lagrange multipliers) X a in the relations expressing l n+1 as a function of Whether all 

these Lagrange multipliers remain arbitrary or not depends on how many of the constraints are 
preserved by time evolution. So we have the same number of constraints as before following 
from the relations 

ye gi n+i \ I 

If these constraints coincide with the ones derived from (|7.12p then we obtain no further 
constraints and the Lagrange multipliers remain arbitrary, as the constraints are automatically 
preserved by time evolution. In general it might however happen that new (secondary) con- 
straints appear or some of the Lagrange multipliers get fixed, see also |18j for an extensive 
discussion. For the case considered here - as the canonical dynamics is just a reformulation of 
the covariant Regge equations - we should obtain the same number of free Lagrange parameters 
as null vectors for the Hessian. This was indeed the case for the examples considered in section 

El 



8 Discrete reparametrization invariant systems 

We so far discussed examples from Regge calculus where for curved solutions we did not find 
an exact form of gauge invariance. One could argue that the reason for this is that in Regge 
calculus one already operates on the space of diffeomorphism invariant geometries, hence one 
would not expect any gauge symmetries to appear. 

If we adopt this view it might be however very difficult to understand the continuum limit. 
In particular for topological theories, such as 3d gravity with cosmological constant it is hard to 
see how in a refinement limit that introduces more and more degrees of freedom (that would all 
be physical) one can obtain a continuum theory, in which only a finite and very limited number 
of degrees of freedom are actually physical. 

Indeed one expects |10t \TT\ [12] that in a continuum limit the (continuum) gauge symmetries 
are restored. The arguments are based on the fact that in a refinement limit the "curvature per 
simplex", that is the deficit angles are getting very small (with the edge lengths getting small 
as compared to the average curvature scale). We have seen that the eigenvalues of the Hessian, 
that correspond to the broken or 'would-be' symmetries go to zero (quadratically) with the 
deficit angles. 

The viewpoint that discretizations lead to a breaking of gauge symmetries that are then re- 
stored in the continuum limit is also strengthened by considering discretizations of parametrized 
(finite dimensional) systems. In those cases the reparametrization invariance is generically bro- 
ken for the discretized systems. One could also argue that one is somehow working on the space 
of gauge invariant data, there is however always a choice of action for which reparametrization 
invariance is fully restored! This shows that it is important to keep in mind that the question 
of gauge symmetry or gauge invariance is one determined by the dynamics - that is the action 
- of the system. 

We will give a short discussion of discretized reparametrization invariant systems, more 
material and applications to numerical integration can be found in [25]. 

To obtain a continuous action with reparametrization invariance we can start from a regular 
Lagrangian L(q, q) where q denote the configuration variables. To obtain a reparametrization 
invariant action we add the time variable t to the configuration variables and use s as an 
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(auxiliary) evolution parameter instead. It is then straightforward to verify, that 

S= H L(q£)t'ds (8.1) 

J Si 

is indeed invariant under reparametrizations s = f(s) of the evolution parameter and the in- 
duced change t(s) = t(f~ 1 (s)),q(s) = g(/ _1 (s)) of the evolution pathes. Here a prime denotes 
differentiation with respect to s. 

One family of discretizations of this action is given by 

N-l 



n=0 

with 

n+1 _ n 

S in ,n + l) ■■= (t n+1 ~ t n ) L((l - a)q n + aq n+l , \ n+l _ Q - ) . (8.3) 

The (Euler-Lagrange) equations of motion are then given by second order difference equations 
which can be obtained from 

q _ ^(n-l,n) ^<%i,n+l) 

° " -^q^ + ^q^- ■ (8 ' 4) 

As an example consider the discretized reparametrization invariant harmonic oscillator with 
a discrete action for one time step given by 

= - \ w (h n + b n+1 ) 2 ) (* n+1 - n (8.5) 

where w parametrizes the strength of the potential term. 

Assume we consider the smallest boundary value problem with prescribed data at time steps 
(n — 1) and (n + 1). A necessary condition for the system to be reparametrization invariant is 
that the Hessian of the action with respect to the variables at the time step n has a null vector 
if evaluated on solutions. Now basically the same discussion as in section [7] will show that if this 
Hessian has a null vector, this is also the case for the matrix 

K := I X dqn+1 X 9tn+1 I • (8-6) 
dt n dq n+1 dt n dt n +' L 



(To consider this matrix has the advantage that we do not need to determine the solution.) 
An explicit calculation gives for the determinant of K 

1 (n n+1 — n n ) 2 1 
detW = A W (U - V? + 4 w2( ^" +1 + ^ • (8 ' 7) 

that is a non-vanishing result. Hence the discrete system does not display (local) gauge sym- 
metries. Note that for a free particle, w = 0, reparametrization invariance persists for the 
discretized theory. As we will see later on the reason is that the dynamic for a free particle 
defined by the discretized theory coincides with the dynamic of the continuum theory. 

This raises the question whether for instance similarly to the case of 3d gravity with cos- 
mological constant discussed in section [6j the gauge symmetry and therefore constraints persist 
for a first order truncation in w. It will turn out that this is not the case. On the other hand 
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one can define the time continuum limit without any problems and recover reparametrization 
invar iance. 

To define time evolution in a canonical formalism we again use the discrete action for one 
time slice as a generating function: 

„ ^(n,n+l) _ Ul q ) + 1 ( , „ + l + _ e) 



dq n ( t n+l _jn) ^ 2^2 

dS {n , n+1) _ l ( g n + l_ g n )2 

<% n ~~ 2(t n+1 -t n ^ 2 - - ! 2 ' 



+1 = ^(n,n+l) = [q q _ i i +1 i +1 _ 

n+1 ^(n,n+l) _ 1 (q n+1 - q n ) 2 1 (l n+1 l n)2 , , 

For the free particle to = the first two equations in (|8.8|) cannot be solved uniquely for 
and one is left with a one-parameter ambiguity. Accordingly there is one constraint 

C w=0 =p? + ^) 2 = (8.9) 

between the canonical data at any time step n. For w > there are at most discrete ambiguities 
in the solutions (q n+1 , t n+1 ) as functions of the canonical data at time step n. To check whether 
we obtain a constraint in a first order truncation in w we expand 

q n+l = O q n+l + w l q n+l ^ ^+1 = Ofl+1 + w 1^+1 

where °t n+1 = t n + ( Pq) (°q n — q n ) is the solution for the w = dynamics. We use these 
expansions in the defintion of the momenta Pq,Pt in (|8.8j) and neglect any higher than first order 
terms in w. From the resulting equations 

„ _ o n _ ( 1 9" +1 -g") 2 ( 1 g n+1 +9")+4(°^) 2 ( ^ H n - V) 

Pq ~ Pi W 4V?(V +1 "f) 

p? = ^ ^ 1 ^ g ?-± (8.11) 

we do not obtain a constraint. Rather these can be solved for instance for the variables l q n+1 and 
°q n+1 (with H n+1 remaining undetermined). That is these (first order) equations fix the values 
for the zeroth order variables °q n+1 , o^n+i were only determined up to a one-parameter 

reparametrization invariance by the zeroth order equations. This behaviour differs from the first 
order dynamics in the cosmological constant for 3d Regge calculus found in section [6l 

On the other hand it is possible to find a continuous time limit, in which reparametrization 
invariance is restored. (This was not possible for 3d Regge calculus with cosmological constant.) 
Again expanding 

q n+l =q n + £ l g n+l ^ = f n + £ (g^ 

we will find that up to second order terms the relation 

C e = P? + \(p n q f + \vj{q n ? + 0{e 2 ) = (8.13) 

between the canonical data. Indeed (|8.13p coincides with the continuum constraint for the 
parametrized harmonic oscillator. 
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Finally we want to point out that for such mechanical systems, in which we have only a 
one-dimensional reparametrization invariance in time direction, it is straightforward to define 
a discrete dynamics which preserves reparametrization invariance. The idea is that the discrete 
system should exactly reproduce the dynamics of the continuous system. To achieve such a 
dynamics one can start with the continuum action and solve for all variables t(s),q(s) except 
for those at some discrete subsets of the evolution parameter s n . Reinserting the solutions in 
the action will result in a sum of Hamilton-Jacobi functions 

N-l N-l rs +1 

S exac t = £4/ s " + V,gV n+ ^ n+1 ) = E / dsC(t(s),q(s)) . (8.14) 

where C is the Lagrangian of the reparametrization invariant continuous system, and t(s),q(s) 
are solutions with boundary data (t n , q n , t n+1 , q n+1 ). 

If we define the momenta for this discretized dynamics according to (18. 8h we actually obtain 
the same definition as for the continuums dynamics 



Pa 



/ 71 71 -4- 1 \ 

dS ( Hj' S ) (t n ,q n ,t n+1 ,q n+1 ) 
dq™ 

j(s n ,s™+ ) / , n n ,n+l n+l\ 



P ? = — - ' ■ (8.15) 

Hence the same constraints between the momenta and configuration variables as for the contin- 
uum dynamics have to hold, which also shows that reparametrization invariance is preserved. 



9 Discussion 

We have shown that we do not have exact gauge symmetries for non-flat vertices in discretized 
4d gravity as defined by the Regge action. In a canonical formulation for Regge calculus this will 
lead to pseudo constraint^) involving data at different time steps instead of proper constraints 
involving data at the same time step. 

This fact should be considered in attempts to connect covariant and canonical quantum 
gravity models, for instance loop quantum gravity, based on proper constraints, and spin foam 
models, that can be seen as quantizations of Regge calculus, which rather lead to pseudo con- 
straints. Indeed the methods discussed in this paper allow to introduce a canonical formalism 
that completely matches the dynamics of the covariant system and moreover using tent moves 
makes a local analysis possible. These techniques seem therefore be ideally suited to be applied 
to other discretized actions, as the discretized Plebanski action [32] used in spin foam models, 
or the area-angle action for Regge calculus introduced in [33] . This could help to derive a closer 
connection between spin foam models and loop quantum gravity. 

In the case that one starts with an action with broken symmetries and obtains pseudo 
constraints, there might nevertheless exist certain limiting cases in which these turn into proper 
constraints. This was the case for the first order dynamics in the cosmological constant for 3d 
Regge calculus. An analysis for 4d Regge calculus for such limiting cases will appear in [28] . 
Starting from these proper constraints it might be possible to extend the constraints, such that 
in the end one obtains a system with an alternative dynamics with exact gauge symmetries. 

The example of reparametrization invariant systems shows that the breaking of the symme- 
tries in discrete theories is due to deviations of the discrete dynamics from the continuum one. 

6 Note that for instance for four-valent tent moves or Cauchy surfaces formed by the boundary of a 4-simplex 
one will obtain proper constraints restricting the physical degrees of freedom to zero, as the dynamics supported 
by these structures is just given by 4d flat space, see [14] . 
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We discussed examples, where different discretizations of the action lead to a symmetry breaking 
and symmetry preserving dynamics. Hence one should be careful in deriving general conclusions 
for quantum gravity models from the fact that the symmetries are broken for the Regge action. 
Alternative actions might display exact gauge invariance, see also |34| for attempts to obtain 
such actions. 

For deparametrization invariant systems we have seen that such an alternative action can 
be always defined by starting from the continuum action and integrating out infinitely many 
degrees of freedom, that is in other words by applying a renormalization group transformation. 
We will explore this idea in more detail in [30]. Indeed in lattice field theories the so-called 
perfect actions, which reflect exactly the continuum dynamics, can be constructed as fixed 
points of the renormalization group flow [35j. These lead typically to non-local couplings (i.e. 
couplings extending over non-adjacent lattice sites). It would be therefore interesting to extend 
the analysis of this paper to non-local actions. 

Note that another property for the exact discretized action (I8.14p is the independence under 
refinements (for the partition function or dynamical observables). Indeed as the dynamics 
already reflects the continuum dynamics, there is nothing gained by refining the lattice. Similarly 
one would expect that a perfect action for discretized gravity - should it exist - would lead to 
triangulation independent partition functions. This is well known for topological field theories 
(such as 3d gravity), where however the actions are still local. For theories with propagating 
degrees of freedom non-local couplings might arise (even at the classical level). 

An alternative way to obtain triangulation independence for instance for spin foam models, 
is to sum over triangulations, which would also include infinitely refined triangulations. This can 
however be only defined at a formal level, for instance by invoking group field theory methods 
[36]. In some sense the construction of a perfect action can be seen as taking the effects of a sum 
over triangulation into account and might in this way circumvent divergencies that appear in 
a sum over triangulations. An investigation into the properties of such actions seems therefore 
promising to us. 

A An example for physically different solutions with the same 
boundary conditions 

In this section we will construct an example of two physically different Regge solutions for the 
same boundary data. The construction uses two tent moves. In the end we will obtain two 
solutions, one flat and one with curvature. As the final triangulation has an inner vertex we will 
check the Hessian of the action evaluated on these solutions and find null eigenvalues for the flat 
solution and only non-vanishing eigenvalues for the curved solution. 

A.l One tent move 

We consider the following four-dimensional triangulation: start with a regular tetrahedron where 
all edges have unit length. Apply a 1-4 move, i.e. add another vertex A, which is connected 
to each of the four former vertices V±, . . . , V4. The length of each of the edges from A to Vi is 
a > 0. This three-dimensional triangulation serves as a starting point for a tent move applied at 
A. The resulting four-dimensional triangulation has yet another vertex B, which is connected 
with each of the Vi, the length of each of the connecting edges being b. Denote the length of the 
tent pole, i.e. the edge between A and B by t. 

The tent pole is the only inner edge of this triangulation, and its Regge equation reads 
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Figure 11: The initial triangulation. The lengths of the edges connecting any Vi to any other Vj is 1, 
the length of the edges connecting A to any of the Vi is a. 



\- dA ABVi 

2^ — ^ — eABy i = °> 



Of, 



(A.l) 



where AABVi * s the area of the triangle between the vertices A, B and Vi, and eABVi ls the deficit 
angle at the corresponding triangle. Due to our symmetric setting, neither of these depend on 
i, and is solved if either of the following equations holds: 



dt 

€ABVi 



t(a 2 + b 2 -t 2 ) = 
- (a 2 - b 2 ) 2 + 2t 2 (a 2 + b 2 ) 

f _ {a 2_ b 2 ) 2 + 2{a 2 + b 2_ 1)t 2_ t A 1 

"" ' ,<U<<OS 2(a 4 + t 2 + (& 2 -t 2 ) 2 -2a 2 (b 2 + t 2 )) 



0. 



(A.2) 



(A.3) 



The equation (|A.2p has t = t± := V a 2 + b 2 as a solution. The only two positive solutions for 
the second equation are 



t = t± := i ^Aa 2 + 4b 2 - 3 ± \/8a^3 \/86 2 - 3 



(A.4) 



The two latter solutions t± can readily be identified as the two flat solutions pointing in 'for- 
ward' and 'backward' direction. (For b > a the solution i_ is forward pointing, whereas i+ is 
backward pointing.) 



A.2 Two tent moves 



We now again start with the three-dimensional triangulation shown in figure lA~Tl but successively 
apply two tent moves to the vertex A. The resulting triangulation has an inner vertex, B, 
connected to the vertices Vi, the length denoted by b. The final spatial triangulation contains 
a vertex C, the lengths of the edges connecting any of the Vi to C being of length c. For 
simplicity, we choose c = a for the moment. Denote the lengths of the tent pole connecting A 
and B, and B and C by si and S2 respectively. Due to symmetry, the Regge equations for s\, S2 
can immediately be solved to be either of 



Sl,2 = t±, t± 



(A.5) 
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with t± := \/a 2 + b 2 and t± given by ()A.4p . Note that there is still one more equation to solve, 
which is the equation for b, the bulk edge connecting Vi and B. The Regge equation for b reads 

dS dA ABVi OAcBVi . 9 A BV t v 3 . . 

~ ~db = —db~ €ABVi + —db~ ecBVi + 3 —ab~ eBViV > = (A ' 6) 

The equation (|A.6[) will have a different form for each of the six choices (|A.5|) for the lengths 
sx,S2j and result in different equations for b. 

We first consider the cases where si,S2 = i±- As we have already seen, this results in the 
deficit angles eABVi = e CBVi to vanish. Since for any i ^ j one has 

db 2y/W=T [ } 

the only way (nonzero b) for (|A,6h to be satisfied is that also CBViVj = 0, so the overall triangu- 
lation is flat. 

It is now not hard to show - in case that either s\ = i+, S2 = t- or s± = t-, S2 = £+ - that 
eBViVj is identically zero for 6 lying in some finite interval (which, after some investigation, can 
be found to be the range of values such that all triangles in the triangulation satisfy the triangle 
inequalities). This gives a one parameter family of flat solutions for the two cases. In the case 
that either s\ = S2 = t + or s\ = S2 = t-, the only positive solution is 

(A.8) 

which is the point where t+ = t~, so this is just a special case of the flat solution. 

Next we consider si = t+ and S2 = t±. Again, one finds that the deficit angle €BViVj has to 
vanish, resulting in a completely flat triangulation. The only positive solution is 

^ =, (A.9) 



V8a 2 -3 

being the solution of t± = t-, hence again a special case of s% = t + and S2 = t-. The case 
s\ = t- and S2 = t± works analogously. 

There is one case left, which is s\ = S2 = t± = Va 2 + b 2 . Since neither of the deficit angles 
has to vanish now, the Regge equation (|A.6|) involves multiple terms, therefore is genuinely a 
transcendental equation, given by 

a (2* - 3 arccos ^ffijff ) + ^fc ~ 4 arccos ^-^S+a^-i)) ) = ° ' 

(A.10) 

It still can be solved numerically. For a = 4, one can plot the derivative —dS/db w.r.t b, see 
Figure E3 

The zero can be evaluated numerically to be b ~ 0.604458. The deficit angles can be computed 
numerically for this case as well, and are all found to be unequal to zero: 

tABVi = e CBVi ~ -0.615334, e BViVj w 1.84412. (A.ll) 

The four-dimensional volume of the triangulation is computed to be 

V = BVabvmv* ~ 0.180461 (A.12) 

Similarly, the three-dimensional volumes of the tetrahedra can be computed and found to be 
positive, as well as all areas of triangles. Therefore, all generalized triangle inequalities are 
satisfied and the solution defines a geometrical triangulation. 
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Figure 12: The derivative —dS/db as a function of b. 



Note that - with the same boundary conditions a = c = 4, there could always be chosen the 
flat solution si = t + ,S2 = t-, e.g. with b = y/3/8, for definiteness. The corresponding volume 
can be computed to be 

V = SVABViVM ~ 0.232924. (A.13) 
which is different from the curved case. 

This gives an explicit example where the the specification of boundary values^, and a fixed 
triangulation (containing an inner vertex) does not result in a unique solution. Apart from the 
flat solution, there is a curved one as well. Also, this curved solution is isolated, i.e. there is no 
continuous symmetry of the solution, which we now confirm by computing the Hessian of the 
action at the solution. 

In our triangulation there are six edges meeting at the vertex B. Therefore, the Hessian is a 

C 




"w 

A 



Figure 13: The 4d triangulation with inner vertex B. 

6x6 matrix the entries being the second derivatives of the action. Since our solution is highly 
symmetric, with b\ = . . . , 64 = b and si = S2, some of the eigenvalues coincide 

Ai w -137.752, A 2 ~ -8.32992, A 3 ~ -8.32393 
A 4 = A 5 = A 6 « 0.339397 . 

One can see that the Hessian possesses no zero eigenvalue, showing that there is no infinitesimal 
symmetry of this solution. It is instructive to compare this solution with the flat, homogenous 

7 In this case we chose a — c = 4 for simplicity, however the same type of curved solution exist also for a ^ c. 
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solution, which is given by 




(A.14) 



In this case there are four null eigenvectors that can be associated to translations of the inner 
vertex in the four directions. 



B Geometric relations in simplices 

Consider a D-dimensional simplex in D-dimensional Euclidean space or in D-dimensional man- 
ifold of constant sectional curvature k / consisting of D + 1 vertices v%, . . . , vd+i- Denote this 
simplex by (123 . . . D+ 1). Any subsimplex is determined by the subset , . . . , Vi n of the vertices 
which span this subsimplex, and will therefore be denoted as (i±i2, • • • , i n )u Similarly, the notion 
(ii*2 • • • , *»») is denoting the subsimplex which consists of all vertices other than i\, . . . , i n . The 
subsimplices in curved space are defined to be the hypersurfaces with zero extrinsic curvature as 
embeddings in the geometry of the higher dimensional simplex. These are in fact also simplices 
of curvature k. An edge (ij) is then just given by the geodesic connecting Vi and Ujjfl 





Figure 15: A spherical triangle in S 2 is a simplex 
of constant curvature 1. Note that the edges are 
Figure 14: The interior dihedral angle in the tetra- simplices in S 1 , which are the great circles, having 
hedron (1234) at the edge (34) = (12). zer0 extrinsic curvature in S 2 . 

Denote the geodesic lengths of the edges (ij) by Uj. Then the (D + 1) x (D + 1) matrix G with 
entries 



G 



ij 



(B.l) 



where the function c K (x) is defined by 

c K (x) := < 



COS (\JTix) k > 

cosh (V — KX ) k < 
is called the Gram matrix of the simplex. For k = we define 



'jl 



D+l. 



(B.2) 



h:,l 



The order of the ik determines also an orientation of the subsimplex, which however plays a minor role in 
this article. 

9 Note that if k > then there are at least two such geodesies. In fact, there are several different curved 
simplices which have the same vertices. For definiteness we will always choose the shorter of the two geodesies, 
and not consider cases in which there are infinitely many such geodesies, which are degenerate anyway. 



28 



to be the affine metric of the simplex, see for instance [23]. For k / we denote by c^j the 
ij-th cofactor of G, i.e. the determinant of the matrix obtained by removing the i-th row and 
j-th column of G. For the flat case n = 0, we take Cjj to be the affine inverse metric, that is 

Then the interior dihedral angle 9^ opposite of the edge (ij) is given by [37J 123! 

cos %i = - 7=n= - ( B - 3 ) 

Denote the volume of the subsimplex spanned by all vertices except Vi and Vj by V^y For 
variations 5 of the geometry of a simplex the Schlafli-identity [38J 

E%)<%1 = P- 1 )«* v (i2...i3+i) ( B - 4 ) 

i<j 

holds. 

With the help of this identity we can derive some formulas needed in section [6J First of all 
if for a given tetrahedron we define the function 

S - = E%)%] " 2kV (12...d+i) (B.5) 

i<j 

we can use it as a generating function for the dihedral angles 



(B.6) 



where Vfe\ is the length of the edge between the two other vertices k, I ^ i,j of the tetrahedron. 
Since partial derivatives commute we can conclude 

which is used in section [6l for instance to see that the constraints derived there are Abelian. 

To compute the derivative of the dihedral angles with respect to the sectional curvature n, 
also needed in section [6l we note from the definition (|B.2|B.3|) of the dihedral angles that 

d6 [kl] _ \- d6 [kl] 1 ,„-! ^ 

_ y d hi\ 1 -i„ 

= ^ < B - 8 > 

ov (ki) 

where in the second line we use (|B.7[) and in the third line the Schlalfi identity (|B.4p . 

The derivative of the dihedral angles can be also computed from their definition (|B.2IB.3|) 

d°[kt\ 111^/ c M c kl \ dG hm 

-7T r = ; — J. , -Z > , CkhCml + CkmChl CkhCkm QfeQm ~ 7T, , {O-V) 

dlop sin 0py] ^CkkCu 2 ^ V c kk Qz J dlop 

see also |23j for simplifications in the Euclidean case. 
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Finally the volume of an Euclidean Z)-simplex can be computed by using a generalization 
of Heron's formula 



VA 



(-1) 



D+l 



1 



det 



1 
1 



1 



r 2 

'12 



I 2 



72 ;2 

to 



'(12.. .D+i) - v -v 2 D (Z)!) 2 

V 1 'lD+1 l 2D+l 

and one can show that for the derivative of this volume we have 

D 2 ^ n D 2 



COS t 



1 

I 2 

1 1D+1 
I 2 

l 2D+l 



^D+1D+1 / 
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(B.10) 
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